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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  of  this  work  was  to  examine  the  frequency  dependencies  of  target  high¬ 
lights  in  order  to  determine  where  computational  savings  could  be  realized  in  broadband 
applications. 

RESULTS 

Significant  savings  can  be  achieved  by  simply  separating  all  frequency  independent  calcu¬ 
lations  from  the  frequency  loop.  In  some  cases,  highly  accurate  approximations  will  provide 
additional  savings.  In  other  cases,  and  contingent  on  the  user’s  fidelity  requirements,  coarser 
approximations  will  yield  yet  more  savings. 

RECOMMENDATIONS 

Separation  of  frequency-independent  factors  and  accurate  approximations  should  be  im¬ 
plemented  as  soon  as  feasible.  Coarser  approximations  should  be  implemented,  but  their  use 
should  be  contingent  on  a  tolerance  test. 
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INTRODUCTION 


Accurate  target  models  that  predict  submarine  and  surface  ship  echo  time  history  are  essential 
tools  for  antisubmarine  warfare  and  ship  vulnerability  studies.  A  realistic  model  must  account  for  the 
target’s  many  highlights  as  well  as  insonification  conditions.  It  must  predict  the  highlight  amplitude 
and  phase  as  well  as  temporal  and  spectral  behavior.  These  predictions  require  that  researchers 
properly  identify  scattering  mechanisms  and  incorporate  them  into  the  model. 

Over  the  past  30  years,  Naval  Command,  Control  and  Ocean  Surveillance  Center  RDT&E 
Division  (NRaD)  researchers  developed  Uniform  Theory  of  Diffraction  (UTD)  target  models  that 
provide  an  invaluable  asset  for  torpedo  signal-processing  algorithms  and  torpedo-simulation 
development.  NRaD  has  exported  these  models  throughout  the  Navy.  Previous  NRaD  target  model 
formulations,  dictated  by  past  weapons’  requirements  and  test  scenarios,  were  monostatic,  far  field, 
and  narrowband.  They  also  used  a  simple  Doppler  implementation.  These  restrictions  simplified  the 
numerical  implementation  (as  well  as  theoretical  development)  considerably.  These  models  serve 
well  in  deep  water  and  other  benign  environments.  However,  the  littoral  environment  requires  more 
sophistication. 

A  promising  avenue  for  study  involves  broadband  processing.  Obviously  simulation  studies  in 
this  area  require  high-fidelity  broadband  target  models.  Practical  considerations  demand  computa¬ 
tional  efficiency.  This  report  examines  the  fi-equency  dependencies  of  target  highlights  and  deter¬ 
mines  how  computational  savings  may  be  realized  without  undo  loss  of  fidelity.  The  report  first 
discusses  general  considerations  followed  by  each  highlight  class,  in  turn. 


GENERAL  CONSIDERATIONS 


Following  Primakoff  and  Keller  (1947),  the  scattering  amplitude  is  taken  to  be  the  product  of  a 
geometrical  factor,  a  phase-cancellation  (directivity)  factor,  and  a  reflection  factor.  In  addition,  we 
include  a  loss  factor  to  account  for  cover  penetration. 


GEOMETRICAL  FACTOR 

For  a  surface,  the  geometrical  factor  is 

G  = 


2r  1  +  —  1  +  — 

i  ^  J 

where  r  is  the  range  and  and  Rg  are  the  principal  radii  of  curvature.  For  a  line, 
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G  = 


ll  + 


2r  sin  y/ 


R. 
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where  y/  is  the  inclination  angle.  This  factor  is  obviously  independent  of  frequency. 


PHASE-CANCELLATION  FACTOR 

The  phase-cancellation  factor  is 

F=P„Py 


where 


±i-r^ 

2 


dr 


is  the  Fresnel  integral  and  the  arguments  are 
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with  k  the  wavenumber.  Now, 

F  ±  (m)  =  C(m)  ±  iS{u) 

where  and  g(w)  are  auxiliary  functions  (Gautschi,  1965).  These  are  odd  functions  and,  for 
0  <  M  <  00 ,  have  the  rational  approximations, 

-/  X  l  +  0.926w 

f(u) « - ^ 

2  +  1.792w  +  3.104w' 

~  2  +  4.1 42w  +  3.492m^  +  6.670u^  ' 
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Alternatively,  series  representations  may  be  used  (Boersma,  1960).  Note  that,  if  ^'min  ~  “/ITmax 
;t.„>4,0.90<|P„|<1.10. 


The  case  of  a  straight  side  must  be  treated  separately.  The  observation  point  will  be  referenced  to 
mid-height.  Let  h  be  the  height  of  the  side  and  y  the  orientation  angle  with  respect  to  the  side.  In 
the  monostatic  case,  it  can  be  shown  (Chabries,  1975)  that  we  may  still  use  our  earlier  definition  of 
Py ,  but  with 


1  , 

rcosy  +—h 


and 


^min 


kh^ 

Alternatively,  the  results  of  Skudrzyk  et  al.  (1973)  may  be  used  when  r  » - .  Here  they  may  be 

'In 

expressed  as 


Pv  = 


-he  ^ 


sm{khcosy^ 

khcosy 


Note  that  the  first  term  in  brackets  is  the  Sommerfeld-MacDonald  factor  for  effecting  the  near  field  to  far 
field  transition  (Ruck  et  al.,  1970).  The  second  term  in  brackets  is  the  directivity  of  a  line  array.  Note 
that  for  khcosy  >30,  Py  is  usually  insignificant. 


REFLECTION  AND  TRANSMISSION  COEFFICIENTS 


Consider  an  acoustic  wave  in  medium  3,  a  fluid  (water),  incident  on  medium  2,  an  elastic  plate  of 
thickness  d ,  which  in  turn  is  backed  by  medium  1,  a  fluid  (water  or  air).  The  reflection  and 
transmission  coefficients  are  given  by  Brekhovskikh  (1980): 


[m'^  -N'^)^  +  \  +  iM  1-^ 


R  = 


ZJ 


1+  ' 

V  -3 


and 


T  =  - 


Z,  p. 
P\ 


z%  V  z^j 
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where 


Z  Z 

M  =  cos^  (26^5  )  cot  P  +  — ^  sin^  )  cot  g , 


jy  _  Z2C  QOS  (2^2^)  ^  (^^2s) 

Z3  sin  P  Z3  sin  Q 


P  =  k2cdcos92c , 


Q=  k2sdcos&2s 


COS0,  ’ 


_  P2^2C 

cos6^c  ' 


7  _  P2^2S 

“  COSP25  ’ 


Z3=^ 

COS  ft 


^2S  ■“  ^ 


2(1 +  <T,) 


The  propagation  angles  (measured  from  the  normal)  are  related  by  Snell’s  law, 

sin ^3  =  kjc  sin^2c  =  ^is  sin  ^25  =  sin^i  . 


kg  N 

Typical  values  for  steel  are  p  =  7.8  x  10^  — j ,  £*  =  2.17  x  10^^  — y  j  ^  -  0-284 .  Hence, 

m  m 

m  m 

=5501  —  and  c.  =3291  — . 

^  s  ^  s 
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Monostatic  Reflection  Factor 

In  the  monostatic  case,  ^3=0  and,  consequently,  all  the  other  propagation  angles  are  zero  as 
well.  Here, 


^2C 


■  +  1  +  / 


V  Z3  y 


V  Z^J 


-^COtP^ 
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with  Pw  =  ky^d .  Then, 
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The  extrema  are  found  from  the  requirement 


d\R 


M\ 


dP, 


=  0  to  be  given  by  the  conditions  cos  P^  =  0  and 


M 


sin  =  0 .  Correspondence  with  maxima  and  minima  depends  on  the  impedances.  Hence, 


R 


A/.extl 


1- 


"2C 


VZ3 


1  + 


VZ3j  Z3 


and 


R 


M,ext2  I 


1-^ 


1.^ 


It  is  useful  to  note  some  impedance  values:  Z^^  =  1.5  x  10® 


=4.1x10’ 

m  s 


4. and 
m  s 


ks 

Zc,steei  =4.3x10’  Therefore,  a  water-backed  plate  has  |PA^„,i„|  =  0  and  \Rj^  ^J^  =  Q.9916, 

while  an  air-backed  plate  has  IP3./ ,„i„|  =  0.6332  and  |P^  „,ax|  =  0.9995 .  For  a  2.5-cm-thick  plate,  the 
minima  and  maxima  are  separated  by  55  kHz. 
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Loss  Factor 


Since  the  penetration  angle  can  range  between  0  and  very  close  to  ~  >  no  simplifications  can  be 
made  to  the  transmission  coefficient. 


HULLS 


Since  the  specular  point  search  and  principal  curvature  calculations  are  independent  of  frequency, 
they  should  be  clearly  separated  fi’om  the  remaining  calculations. 


SAILS  AND  CONTROL  PLANES 


Since  the  specular  point  search  and  principal  curvature  calculations  are  independent  of  fi-equency, 
they  should  be  clearly  separated  from  the  remaining  calculations. 

EDGES  AND  WEDGES 


We  will  use  cylindrical  coordinates  .  The  wedge  surfaces  will  be  defined  as  ^  =  0  and 

(p  =  l7t-  ^wedgc  •  Let  V  be  defined  as 

v:;r  =  2;r  -  ( 


‘'Wedge  • 


Note  that  an  edge  has  v  =  2 .  For  plane  wave  scattering  from  a  straight  wedge,  Pauli  (1938)  gives  the 
monostatic  solution: 

‘I^Paul,  =VsW  +  Va(0), 


where 


>W  = 


2  .  7t 
— sm— 


24^Vv  v\ 


2  cos 


¥ 


_ — 2 - 


¥ 

COS - cos  — 

V  V 


and 

F^{z)  =  ^  e-'^Vr 


IS 


the  complementary  Fresnel  integral.  of  Abraraowitz  and  Stegun  (1965). 


For  a  straight  wedge  of  length,  h ,  and  scattering  of  a  spherical  wave,  the  solution  is 
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where  y  is  the  angle  with  respect  to  the  edge  and  Py  was  discussed  in  the  Phase-Cancellation  Factor 
section. 


HULL  INTERNALS 


The  problem  here  is  in  the  determination  of  the  transmission  coefficients.  Since  the  penetration 
angle  calculations  are  independent  of  frequency,  they  should  be  clearly  separated  from  the  remaining 
calculations. 

SAIL  INTERNALS 

The  problem  here  is  in  the  determination  of  the  transmission  coefficients.  Since  the  penetration 
angle  calculations  are  independent  of  frequency,  they  should  be  clearly  separated  from  the  remaining 
calculations. 


PROPELLERS 


The  hydrofoil  formulation  of  Lengua  (1991)  is  used.  Since  the  specular  point  search  and  principal 
curvature  calculations  are  independent  of  frequency,  they  should  be  clearly  separated  from  the 
remaining  calculations. 


PROPELLER  CAPS 


Consider  a  plane  wave,  of  wavenumber,  k ,  incident  upon  a  finite  cone  with  a  flat  base.  Denote 
the  half-angle  of  the  cone  by  y  and  the  radius  of  the  base  by  a .  The  height  is  then  h  =  a  coty  . 

Polar  coordinates  (r,  (f)  are  used,  with  the  origin  at  the  tip  of  the  cone.  The  values  ^  =  0  and  <j)  =  Tt 
represent  axial  incidence  upon  the  tip  and  base,  respectively.  For  the  purposes  of  this  discussion,  the 

71  71 

base  is  always  shadowed.  Thus,  only  the  interval,  Q<<p<  —  ,  will  be  considered.  When  ^  / 

(broadside),  there  is  specular  reflection  from  the  surface  of  the  cone. 

It  is  useful  to  examine  the  problem  from  a  geometrical  theory  of  diffraction  perspective.  First 
note  that,  for  ()<<f><y  ,  two  singly  diffracted  rays  from  the  edge  are  returned  to  the  source.  Hence, 
the  backscattered  field  is  an  oscillatory  function  of  ka ,  due  to  interference  between  these  rays.  For 

K 

y  <^  <  —  ,  only  one  singly  diffracted  ray  from  the  edge  is  returned  to  the  source.  Hence,  the  field 
does  not  oscillate  as  a  function  of  ka  .  Also  note  that  for  ^  =  0 ,  the  edge  diffracted  rays  have  an 


axial  caustic. 

For  a  unit  amplitude  incident  field,  the  scattered  field  is 

ikr 

where  D  is  the  diffraction  coefficient.  Kouyoumjian  (1977)  has  determined  uniform  solutions  through 
the  axis  and  broadside.  Let 


n  = 


2  7t 


Then,  for  0<(p<Y^ 

D{^)  -  ka - 


.  n 

sin— 


n 


n  ^  11  3n: 

cos — 1  cos - cos — 

n  n  n  , 


where  is  the  zero-order  Bessel  function  and  the  uppier  (lower)  sign  applies  to  the  hard  (soft)  case.  For 


1  ^  2/Mrcscxcos(/+^)+—  /  \  r 

~  ; — e  ^  ^  tan(^  +  ^)r  2ka  esc  y  cos(x 

4  Y^rsin^  *■ 


+ 


where 


r(x)  = 


1  F(x) 

l--e^  forx>0 

2  \  X  2ix 

^  +  _iLiZ  for;c<0 
2  V  X  2ix 


with 


Now,  for  z  «l,  F(z^  »  e 


F{z)  =  2i^‘  dr  . 

sr  ,-,42 


^[m -2ze  * — z^e  ^+--- 
3 


and,  consequently,  for  x  «  1 , 


r(x)  a -i— X  .  Thus,  at  broadside. 


<!-')■ 


(kaf 


.Sk 


9;rsin^  cos;' 
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11  3  1 

For  \z\  »  1 ,  F(z)  »  1  +  —  / - 2-+---  and,  consequently,  for  |x|  »  1 , 

^  z  4  2 

j  /^_£'j  2  2 

rfx)«l - - .  This  is  typically  a  good  approximation  to  within  a  few  degrees  of 

^  2  \x  2  X 

broadside. 

Note  that  there  is  a  discontinuity  at  <!>  =  /■  The  transition  about  this  region  is  extremely  difficult 
to  formulate.  However,  the  scattering  amplitude  there  is  typically  so  small  as  to  not  be  a  concern. 


CORNERS 


Let  us  consider  a  dihedral  comer,  represented  as  a  flange  on  a  cylinder.  It  is  useful  to  first  discuss 
the  special  case  of  a  plane  wave,  of  wavenumber,  k ,  and  unit  amplitude,  incident  upon  a  finite 
cylinder,  of  radius,  a ,  and  length,  2d ,  at  whose  center  is  a  flange,  of  height,  h ,  and  at  a  right  angle 
to  the  cylinder.  Let  ^  denote  the  angle  of  incidence  with  respect  to  the  cylinder  normal.  The 
scattering  amplitude  depends  on  the  extent  of  the  wave  reflected  from  the  comer  back  to  the  source, 
or  “aperture”  of  the  comer.  Now,  monostatic  reflection  from  the  cylinder  and  flange  is  equivalent  to 
bistatic  reflection  from  the  cylinder  to  the  virtual  image  of  the  source.  The  flange  acts  as  the  plane  of 
reflection.  If  the  cylinder  were  of  infinite  length,  the  aperture  would  be  L  =  2h\mi9 .  For  a  finite 
cylinder,  the  aperture  is  limited  to  its  length.  Thus, 

L  =  2miii[c/,/itan^]. 


Bistatic  reflection  from  a  cylinder  of  length,  L ,  has  the  amplitude. 


P  = 

^  Basic 


Jka  cos  9 
\n 


L, 


assuming  kacos9  »  1 . 

For  a  spherical  wave,  from  a  source  at  a  range,  r»2d,  the  backscattering  amplitude  is 


Let  us  now  consider  the  general  case.  Let  the  comer  have  an  angle  ~  ^  and  an  aperture 


covering  the  interval,  x,  <  x  <  X2 .  Then,  it  can  be  shown  (Sides,  1976)  that  the  backscattering 
amplitude  is  multiplied  by  a  correction  factor. 


,2  N 


w=0 
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where 


^  ^cos(^  +  /?„)  sin[  Mx  sin(<^  +  A  )] 


r„  =-2A:(r„-r) 


and 


Pn 


-  arctan 


—  COS(p 

^ _ 


1 

2 

h 

r^cos^ 

1 

U„  J 

=r^  +xl  +2rx„sin^. 


with  x„  =x^+  nAx  and  N  =  —  .  A  typical  increment  is  Ax  =  — . 

Ax  2k 


In  the  continuum  limit, 


B 


t: 


cos  ^  +  arctad 


xcos^  ; 
+  xsin^; 


+x^  +2rxsinj!J 


-2ik\  ^r^^x^+2rxsin^-r] 

€  ^  ^dx , 


which  cannot  be  evaluated  analytically,  the  summation  expression  is  more  suited  for  numerical 
calculation.  However,  if  r  — >  oo, 

sin(H.  sin  -2<*ii|i2.sin  ^ 


ycos^ 


kL  sirup 


Note  that  this  is  the  product  of  three  terms.  The  first  may  be  considered  an  obliquity  factor.  The  second 
is  the  directivity  of  a  line  array.  Hence,  a  non-right  comer  is  analogous  to  non-normal  incidence.  The 
third  term  represents  a  shift  of  the  phase  center. 

When  r  »  |x,_2 1 , 


B 


Ttr  cosp 
2kL^ 


J—irsinp  +  Xi) 

.I^(rsmp  +  Xj)  ] 

1 

\  nr  Jj 

JkrsiTT  ^ 


This  is  usually  a  satisfactory  approximation. 
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TANKS 


Tanks  can  be  very  complex  physical  structures  and  very  difficult  to  model  in  detail.  They  do  not 
wrap  completely  around  the  hull  but  are  separated  by  longitudinal  plates.  They  may  be  reinforced 
with  ribs  and/or  stringers.  Some  may  contain  air  flasks.  They  contain  fluids  of  various  types  and 
may  hold  a  combination  of  sea  water  and  air. 

To  allow  analysis,  a  number  of  simplifying  assumptions  are  made.  The  separator  plates,  stringers, 
and  flasks  are  not  modeled.  Therefore,  only  dihedral  comers  are  considered.  Each  tank  is  uniformly 
filled  with  a  specified  fluid.  A  tank  is  taken  to  be  formed  by  concentric  hulls  bounded  by  deep 
frames. 

Multiple  reflections  (bounces)  within  a  tank  are  designated  by  a  pair  of  integers,  [A/,  Af] , 

following  the  notation  of  Sides  (1976).  N  refers  to  the  number  of  times  a  ray  is  reflected  from 
bulkheads.  M  refers  to  the  number  of  times  a  ray  is  reflected  by  hulls. 

For  the  case  of  plane-wave  incidence  and  right-angle  comers,  the  backscattering  amplitude  of  an 
[A^,  M]  bounce  is 


^Basic 


_  nA^-I  pN  nM-1  pM 
—  TVgj  iVg2  iVoH  - 


kacosO 

An 


]  • 


This  is  composed  of  several  factors.  The  square-root  term  is  the  amplitude  of  an  ideal  comer  reflector,  as 
discussed  in  the  previous  section.  There  is  a  reflection  coefficient,  R ,  for  each  wall  of  the  tank.  The 
factor  in  brackets  consists  of  “transmission  windows”  of  length,  d-^  and  »  multiplied  by  the 
transmission  factors,  7]  and  7^ .  Since  rays  must  pass  through  either  the  leading  bulkhead,  outer  hull,  or 
both,  the  appropriate  factors  must  be  used.  In  terms  of  the  transmission  coefficients,  the  three 
possibilities  are  7^,,  Tq^i,  or  Therefore,  d  must  be  divided  into  regions  of  constant 

transmission  factor. 

The  general  case  includes  a  correction  factor  as  discussed  in  the  previous  section  on  comers. 

One  possible  way  to  account  for  the  effects  of  clutter  is  through  a  mean  free  path  approach.  This 
is  obviously  a  stochastic  analysis,  so  agreement  is  only  expected  in  the  aggregate  and  not  for 
particular  realizations.  The  mean  free  path,  A  ,  may  be  inferred  from  schematics  or  experimental 

-V 

data.  The  amplitude  of  a  particular  bounce  is  then  reduced  by  a  factor,  e  . 

It  is  readily  apparent  that  because  of  the  reflection  and  transmission  coefficients,  the  frequency 
dependence  is  a  complicated  function.  Experience  shows  that  higher  order  bounce  numbers  do  not 
make  an  appreciable  contribution,  so  that  the  tank  response  may  be  simply  taken  as  the  [l,l]  bounce. 
At  this  level  of  approximation,  7fp„  may  be  taken  as  unity. 


FRAMES 


The  problem  of  scattering  from  hull  stiffeners  is  very  complex  and  has  been  the  subject  of 
considerable  debate.  We  use  the  Hayek  and  Karali  (1993)  and  Hayek  (1997)  formulation.  Consider 
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a  plane  wave  of  amplitude,  <I>o ,  incident  at  an  angle,  6^ ,  with  respect  to  the  normal,  on  an  elastic 
plate  of  thickness,  h .  Let  p,  D,  and  G  be  the  density,  flexural  rigidity,  and  shear  modulus  of  the 

2  ^  / 

plate,  respectively,  and  let  k  =  ^ /■^2  • 


D  = 


Eh^ 

12(1 -v/^) 


G  = 


E 

2(1  + v)’ 


with  E  and  v  the  Young’s  modulus  and  Poisson’s  ratio  of  the  plate. 

The  reflection  coefficient  of  an  infinite  plate  without  discontinuities  and  with  fluid  loading  on  one 
side  only  is  given  by 


where 


Here, 


a(^9g)  =  Fj  +  Fj  sin^  ^0 

6(^0)  =  Jsin^  0Q  -  F,  sin^  Oq  +  Fj  jcos^o . 


pco 


1  /l' 

■  +  • 


k^G  \2D 


24  1 

1-  V 


■  +  1 


^2  = 


phoa^ 

ph^co^  1 

_  24  1 

f"'! 

-A 

1 

k^ 

Kq 

Uk^GD  D 

;r^  1-  V 

V  ^0  / 

1 

'  kl  K^Gh 


4V3po 
„2  ^ 


Q 


Pq  ^ 
p  Q 

and  we  also  define 

*  kl  D  2V3  Co  p  ' 


Poty^ 

1  ph^co^ 

4V3  1 

D  Uk^GD_ 

1-V 

UoJ 

'  2V3  Co  Q' 

12 


The  thin-plate  speed, 


and  Q  = 

/(^c 


where  the  critical  frequency  is 


CO 


C 


^/l2co^ 

hCp 


The  result  for  the  scattered  field  is,  after  a  steepest  descent  approximation, 
exp 


i\  Kr-^ 
4 


^>~a>h-V(R  'll— L^— llicos^ _ Co+C,sin^+Qsin  6 _ 

Fj  +  Fj  sin^  ^  +  zcos^sin‘' F,  sin^  Fj  j 


±2;ri^  residues  [between  S  and  L], 


where  6  is  the  observation  angle  and  the  residues  in  this  expression  are  from  the  poles  enclosed  between 
the  steepest  descent  path  and  the  original  path  of  integration.  The  ±  sign  is  due  to  the  fact  that  the  two 
paths  cross  one  another.  Hence,  the  sign  depends  on  where  the  pole  is  located,  that  is,  whether  it  is 
encircled  clockwise  or  counterclockwise.  The  residues  may  be  calculated  as 


residue  of  Uj  =  <l)o[l-F(^5)jexp[zA:orcos(a-0)]cosa 


Q  +  C,  sin  a  +  Cj  sin^  a 


—  Qisina) 
da  ^  ’ 


where  ^sincty  )  =  0  >  with 

Q{X)  =  F,+  - FjA'  +  Fj] , 


and 


da 


g(sin  a)  =  2F3  sin  a  cos  a  +  z  sin  a  -  F2  +  F,  (sin^  a -2  cos^  a^  -  sin"  a +  4  sin^  a  cos^  a  j . 


The  result  is  only  valid  for  large  k^r  ,  and  it  fails  near  grazing  because  cos^  ->  0 .  The  coefficients, 
C„ ,  are  given  by 

■^^iZ jk\F^ 


Imcop^-Z^klFJ^  -Zfk^F^I^ 


■cos  Or, 
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where 


C,= 


427ciZ„ktF^  .  ^ 

- "  °  - sin ^0 cos Oq 

Imcop^  +  Z„k^  FJ^ 


Q  = 


yf^iZj^k^F.^ 


ImcopQ  -Zj-klFJ^  -Zj^k^F^I^ 


cos  6., 


•>-  Q{X) 


Zj  and  are  the  line  force  and  moment  impedances,  respectively,  of  the  frame. 

It  is  obvious  that  this  has  an  extremely  complicated  frequency  dependence  before  even 
considering  the  dependencies  of  Z^  and  .  There  is  no  satisfactory  way  to  approximate  these 

results  over  any  appreciable  frequency  range.  However,  it  is  possible  to  achieve  significant 
computational  savings. 


APPROXIMATIONS  TO  INTEGRALS 

Calculation  of  the  /„  via  numerical  integration  involves  considerable  computation  time.  For  high 
frequencies  they  may  be  approximated.  The  procedure  is  as  follows.  Write 

Q{A)  =  -  Zi'lz'  -  Zj']  +  F5  +  F3Z' , 

where 

V=i{F,±VF,"-4/;}. 


Note  that  for  Q  »  1 ,  the  and  terms  are  small.  Also  note  that  for  most  materials,  F^  »  F2  (and 

F3  »F,).  Hence,  the  integral  is  dominated  by  the  regions  about  =  Z2^ .  The  approximation 
consists  of  using  this  value  of  Z  within  the  radical,  so  that 


/.»f 

Jl-ao 


Z"JZ 


A*  - 


f  \ 

F-  +F  -i  — 

1' 

V 

Observe  that  for  odd  n,  I„  =  0 .  Let  Z  =  -Jx  ,  then  for  even  n , 


i  _i 
^dx 

{x-z,){x-z,)’ 
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where 


^..2  = 


F,+/ 


I 


F+/  ^ 


F2-i . 


\ 


Now, 


f; 


^dx 


for  |arg;ff|  <  tt ,  [arg/l  <  ;r,  and  0  <  Re//  <  2  (Gradshteyn,  1980).  Thus, 


TT 


(-z,)j'->-(-z,)i 


<-0 


CSC 


(22+1) 


/r 


Specifically,  we  have 


7o«- 


1  1 


and 


4 


Let  us  consider  a  5-cm-thick  steel  plate  (coincidence  frequency  of  4.7  kHz).  The  “exact”  and 
approximate  results  are  compared  in  figures  1  and  2.  As  can  be  seen,  the  general  behavior  is 
followed,  but  not  the  oscillations.  This  is  not  as  serious  a  problem  as  might  first  be  expected, 
because  their  effect  is  mitigated  by  the  manner  in  which  they  enter  the  C„ .  Figures  3  through  1 1 
compare  the  results  for  the  non-angular  parts  of  Q ,  C, ,  and  Cj  for  three  cases.  Case  1  has 

Zf  =10U  and  Z„  =  lO^  N-s.  Case  2  has  Zf  =10®/  and  =  10^/N-s.  Case  3  has 
m  m 


Z.  =10'/^ 
m 


and  Z„  =  10®/N-s.  Some  general  observations  may  be  made,  based  on  cases  ran 


between  the  ones  presented.  Q  shows  good  agreement  for  Q  >  1  in  all  cases.  C,  shows  the  general 


behavior,  but  not  the  oscillations.  Errors  will  be  less  than  3  dB  for  Q  >  2  ,  except  when  |Z„|  >  10® . 
Cj  shows  good  agreement  for  Q  >  6  in  all  cases.  It  shows  good  agreement  for  Q  >  2  ,  except  when 
|Z^|>10®. 
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300 


Figure  1.  Comparison  for  1^. 


Figure  2.  Comparison  for  . 
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lo  ol 


Figure  3.  Comparison  for  non-angular  part  of  Q :  Case  1 


Z_f  =  10''3  i  and  Z_m  =  10''3  i 


Figure  4.  Comparison  for  non-angular  part  of  C, :  Case  1 
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Figure  6.  Comparison  for  non-angular  part  of  :  Case  2 


0.34 


Figure  7.  Comparison  for  non-angular  part  of  Q :  Case  2. 


Z_f  =  10'^6  i  and  Z_m  =  10'^3  i 


Figure  8.  Comparison  for  non-angular  part  of  Cj :  Case  2 


U  ol 


Z  f=10'^3ianclZ  m  =  10^i 


2  4  6  8  10  12  14  16 

Normalized  Frequency 


Figure  9.  Comparison  for  non-angular  part  of  Q :  Case  3. 


Z_f=10''3iandZ_m  =  10'^i 
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Z_f  =  10^3  i  and  Z_m  =  lO'^  i 


Figure  1 1 .  Comparison  for  non-angular  part  of  C2 :  Case  3. 


REGIONS  THAT  MAY  BE  APPROXIMATED 


Note  that  the  scattering  amplitude  is  significant  only  in  a  small  region  about  the  coincidence 
angle,  6^ ,  where 


sin^^ 


1 

Vq 


The  width  of  this  region  decreases  with  increasing  frequency.  For  most  purposes,  the  amplitude  can  be 
treated  as  zero  when  more  than  6  degrees  from  0^ .  In  any  case,  the  integral  approximations  are  more 
than  sufficient.  Thus,  numerical  integration  is  only  required  in  the  immediate  vicinity  of  6^ . 


SUMMARY 


The  frequency  dependencies  of  target  highlights  have  been  examined.  There  are  many  areas 
where  computational  savings  may  be  obtained,  without  loss  of  fidelity.  Approximations,  within 
tolerable  errors,  leading  to  further  computational  savings  have  been  identified.  These  should  allow  a 
practical  broadband  simulation  to  be  realized. 
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